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Abstract 

We propose a reduced form set of two coupled continuous time equations linking the price 
of a representative asset and the price of a bond, the later quantifying the cost of borrowing. 
The feedbacks between asset prices and bonds are mediated by the dependence of their 
“fundamental values” on past asset prices and bond themselves. The obtained nonlinear 
self-referencing price dynamics can induce, in a completely objective deterministic way, the 
appearance of periodically exploding bubbles ending in crashes. Technically, the periodically 
explosive bubbles arise due to the proximity of two types of bifurcations as a function of 
the two key control parameters b and g, which represent, respectively, the sensitivity of the 
fundamental asset price on past asset and bond prices and of the fundamental bond price 
on past asset prices. One is a Hopf bifurcation, when a stable focus transforms into an 
unstable focus and a limit cycle appears. The other is a rather unusual bifurcation, when 
a stable node and a saddle merge together and disappear, while an unstable focus survives 
and a limit cycle develops. The lines, where the periodic bubbles arise, are analogous to 
the critical lines of phase transitions in statistical physics. The amplitude of bubbles and 
waiting times between them respectively diverge with the critical exponents 7 = 1 and 
v = 1 / 2 , as the critical lines are approached. 
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1 Introduction 


Much has changed since the last 10 years, when the so-called “great moderation” was regarded as 
the biggest success of modern economic theory, and research had found that the term “bubble” 
was not an attractive explanation for the lack of quantitative understanding of real estate prices 
PQ. Now we live in the “great experiment” era of central banks in which a plethora of “code red” 
policies [2], such as quantitative easing, zero interest rates, large-scale asset purchases (leading 
to de facto currency “wars” and their potential debasement) are being put forward, one after 
another, without any certainty about the social and economic consequences that these policies 
might bring about. The “great experiment” will likely continue for one, perhaps two or even 
more decades (in different forms) as in the Japanese scenario since 1990 and the bursts of their 
great market and real-estate bubbles. Many new bubbles in part fueled by these policies are 
likely to emerge. 

In this context, the present article investigates the impact of feedback loops between asset 
prices and bond prices (or between monetary base and interest rate spreads) that lead to regimes 
of periodically explosive bubbles followed by crashes. While asset and bond prices tend to mean 
revert to their respective fundamental prices, we propose a simple framework to account for the 
fact that the latter are not fixed but depend themselves on the state of the economy represented 
by the instantaneous asset and bond prices. This results in a rich nonlinear coupled dynamics of 
asset and bond prices that we explore. Our work is related to a rich literature on bubbles and 
crashes, which has documented that, surprisingly from the point of view of the efficient market 
hypothesis, they appear repeatedly with remarkable regularity [3j. 

There exists a number of models advanced to rationalize the appearance of these events. 
A common explanation to these episodes is to regard them, following Keynes [4], as outbursts 
of irrational speculations , such as manias, panics, and similar herding and crowd behavior 
Such a collective behavior can be interpreted as leading to phase transitions, the bubbles bursting 
in crashes, in an ensemble of interacting agents mu- 

Another approach is based on multiple-equilibria models [T2l IT5] . where several equilibria 
arise due to the existence of different types of traders, including dynamic hedgers. Then random 
motion can lead to jumps between different equilibria, which imitates crashes. Similar jumps can 
occur in markets with inexperienced traders [16]. In these models, in the absence of dynamic 
hedgers or inexperienced traders, no bubbles and crashes occur. The nonlinear dynamical model 
of a market with speculative traders Era can also exhibit several equilibria, with jumps between 
different regimes caused by random noise. 

In liquidity shortage models [18], a crash can occur when price plummets due to a temporary 
occasional reduction in liquidity. Bursting bubble jnodels assume a scenario when all market 
traders realize that an asset price is larger than its fundamental value, but they keep buying the 
asset in the belief there are others who do not yet know that the asset is overpriced, and to whom 
they expect to sell the asset at a higher price. At some random point of time, everyone (or a 
sufficient large fraction of investors) realizes that too many are aware of the overpricing, which 
results in the bubble burst and a crash [T9U23] . 

In lumpy information aggregation models, the existence of overpricing is known only to a 
fraction of the traders, while others are unsuspecting. At some point, the uninformed traders 
suddenly discover the existence of overpricing, which leads to a sharp decline in prices imitating 
a crash [241 - f26j . A sudden change in information can also produce a jump in prices [2?] . 

These different models thus explain booms and crashes as caused either by irrational specu¬ 
lations, or by sudden variations of some market characteristics, such as asymmetric or subjective 
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information. 

In the present article, we abstract from the specific micro-mechanism to emphasize the role of 
delays and of nonlinear feedbacks [28] . We propose a reduced form set of two coupled continuous 
time equations linking the price of a representative asset and the price of a bond, the latter 
quantifying the cost of borrowing. We argue that the presence of simple feedbacks between 
asset prices and cost of borrowing can induce, in a completely objective deterministic way, the 
appearance of periodically exploding and collapsing bubbles ending in crashes. This mechanism 
does not require heterogeneous agents, or sudden changes in system characteristics, or multiple 
equilibria. 

The plan of the article is as follows. In Sec. 2, we formulate the dynamical model characteriz¬ 
ing the relation between an asset price and a bond price. In Sec. 3, we give the detailed analysis 
for the existence of stationary solutions and their bifurcations. The various temporal behaviors, 
exhibiting bubbles and crashes, are studied in Sec. 4. And Sec. 5 concludes. 

2 Model formulation 

We consider an economy represented by one representative asset and one bond quantifying the 
cost of financing. Our goal is to derive simple dynamical equations governing the interplay be¬ 
tween asset and bond prices, which emphasize the existence of highly nonlinear unstable regimes 
associated with booms and crises. The interdependence between stocks and bonds is one of the 
most important factors to asset allocation decisions. It also reflects the impact of central banks 
monetary policies on the growth of economies. For instance, Knut Wicksell pointed out that the 
greater the difference between the ‘natural’ interest rate and the ‘market’ rate, the bigger the 
subsequent booms and bursts [29]. 

A unified framework for pricing consistently all assets is expressed in terms of a stochastic 
discount factor (SDF) M(t), also called pricing kernel, pricing operator, or state price density [30]. 
In this framework, the price of an asset is equal to the sum of expected future payoffs discounted 
with the SDF, which embodies all the macro-economic risks underlying each security’s value. 
Under an adequate definition of the space of admissible trading strategies, the no-arbitrage 
condition translates into the condition that the product of the SDF with the asset price p{t) of 
any admissible self-financing trading strategy, implemented by trading on a financial asset, must 
be a martingale, 

M(t)p(t) = E [M[t')p{t')} , (1) 

where t' > t refers to a future date and t is the present time. The expectation E[.] is taken with 
respect to all the available information up to time t. Introducing the discount factor from t to t' 
defined by m(t,t') := expression (HD becomes 

p(t) =E[m(t,t')p(t')] . (2) 

The discount factor m(t, t') also characterizes the cost of financing [3Tj . 

Let us consider the zero-coupon bond B(t, t + r) that matures at time r in the future. We fix 
r to represent the typical investment horizon of a representing investor, say three or six months. 
By definition, B(t + r, t + r) = 1 and, thus, because of Eq. ((2D, one has 

B(t,t+ r)—E[m(t,t+ r)] . (3) 

The bond price B(t,t + r) thus embodies the information on the discount factor and the price 
of financing in the economy. 
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We assume that there is a natural equilibrium price p* for the asset and B* for the bond 
value considered above. Close to this equilibrium, we assume that p{t) and B(t, t + r) obey the 
simplest two-side exponential mean reversal dynamics towards the equilibrium: 


dp p 2 dB B 2 

dt ^ p* ’ dt B* 


(4) 


where we call B(t) := B(t, t + r) for simplicity of notations. The constant coefficients in the right- 
hand-side of these two equations have been absorbed in a rescaling of time and in the definition of 
p* and B*. These equations are known as logistic equations, and enjoy the following properties: 
(i) p(t) and B(t) converge, respectively, to p* and B*\ (ii) the convergence is exponentially fast, 
whatever the initial values, smaller or larger than the fixed points p* and B*. This represents a 
convenient concise description of a mean-reversal dynamics. 

In reality, p* and B* are not fixed but depend themselves on the specific state of the economy, 
which we can characterize by the instantaneous value p(t) and B(t). Intuitively, as the bond price 
increases (interest rates decrease), the price level tends to grow due to cheap financing. As the 
price increases, this tends to push later the rate upward as access to funding competes with 
investments in the appreciating asset and as central banks attempt to cool down a possible 
over-heating of the economy. As the price accelerates, the bond price may decrease drastically, 
the interest rates shoot up, and then a price crash may ensue. To capture this intuition, and 
given the structure of expressions (J 2 J) and we assume that p* is a function of the product 
7T (t) := p(t)B(t), while B* is just a function of p(t). 

The explicit dependence of p* on n(t) and of B* on p(t) can be derived in two ways producing 
the same result. The simplest way is to assume that the relative variation of p* is proportional 
to the variation of n(t): 

Sp* oc p*Sn(t) . (5) 

Similarly, the relative variation of B* is proportional to the variation of p(t): 


SB* oc B*Sp(t ) . 


( 6 ) 


This yields 


p*(n) = p*{0)e bn B*(p) = B*(0)e 9P , 

where the parameter b is a fundamental log-price rate and g is a log-discount rate: 

d In p* d In B* 


b = 


r>Tr 


9 = 




(7) 

( 8 ) 


Without loss of generality, it is possible to reduce the evolution equations (J3| with (J7J) to dimen¬ 
sionless normal forms pITl CT4| . Denoting 


x(t):=p(t), z(t) := B(t) , 


(9) 


we obtain the equations 


dx 2 

— = x — x e 
dt 


—bxz 


dz 

dt 


z 2 e~ 9X 


( 10 ) 


In what follows, we keep in mind this reduction and that all quantities are dimensionless. By 
definition, the asset and bond prices are non-negative x > 0 and z > 0 . 
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The structure of the coupled nonlinear logistic equations fllOp is reminiscent of the dynamics 
of symbiotic biological species with population fractions x and z, influencing each other through 
their carrying capacities [32H3^j . or of the predator-prey dynamics with mutual interactions 
between predators and preys [35]. Another analogy is the co-evolution of different groups in 
a structured society [36|l37]. However, the basic point in our evolution equations is that the 
functions x and z are not different species, but the asset and bond prices in a market. 

When interest rates decrease, the effective bond price increases, increasing the effective price 
xz, and the fundamental price level p* grows. Then, since p* ~ e bxz , this implies that the rate b has 
to be positive or, more generally, non-negative to represent a realistic economic regime. In turn, 
the increase of the asset price x tends to trigger eventually an increase in the cost of investments, 
which tends to push down the effective bond prices, which means that the coefficient g is usually 
negative (non-positive). Thus, the standard situation in a financial market corresponds to the 
conditions 

b > 0 , g < 0 . (11) 

The other way of deriving the equations fllOH is to assume that n(t) and p(t) are initially 
small so that p* and B* can be represented in the form of Taylor expansions 


P = 


E 

n =0 


TC 


n\ 


b " = E 


n =0 


-j p n 
n\ 


( 12 ) 


Then, it is possible to extrapolate the above expansion to finite values of n(t) and p(t) by invoking 
resummation techniques, for instance self-similar approximation theory [38TH(Tj . In the course of 
the resummation, it is necessary to impose the restriction of non-negativity of p*(t) and B*(t) |4Ij. 
As a result of such a self-similar resummation, under the condition of semi-positivity of p* and 
B*, we come to the same form (J7J) corresponding to a self-similar exponential approximant ||42j]. 

The above derivation of the coupled dynamics (flOl) of asset and bond prices does not include 
any stochastic component as we consider averaged macro-quantities in order to emphasize the 
nonlinear feedback loops between asset and bond prices. In reality, a stochastic structure would 
need to be added, which can be both additive and multiplicative (corresponding to parametric 
noise). Here, we focus on the deterministic structure of the dynamics to unravel the main 
consequences of the feedback loops. In particular, we uncover a regime of periodically collapsing 
bubbles with remarkable properties. 


3 Stationary solutions of evolution equations 


Equations (HU possess the following trivial fixed points: {0,0}, {1,0}, and {0,1}. In addition, 
there are nontrivial fixed points, which are the solutions of the equations 

x* = exp (bx*z*) , z* = exp(gx*) . (13) 


The stability of each fixed point is described by the characteristic exponents 


A 



bx* z* 


2 ± x*\/bz*{Ag + bz*) 


(14) 


A fixed point is stable (respectively, unstable) for negative (respectively, positive) real parts of 
the exponents. The trivial fixed points {0, 0}, {1, 0}, and {0,1} are always unstable for all values 
of b and g. 


5 





Depending on the parameters b and g, np to three nontrivial fixed points {x{,z*}, {x^z^} 
and {£ 3 , 23 } can exist, which we enumerate by descending values x\> x^> x 3 of the reduced 
asset price variable. According to condition ([IT]) , we consider b > 0 . Although the parameter g 
should be non-positive, we also study the regime with g positive in order to exemplify the analytic 
continuation of the fixed points through the boundary g — 0. This analytic continuation allows 
us to explain the dynamic behavior of the solutions when approaching this boundary g — 0 . 

The b-g plane is partitioned into four regions by the lines g — 0, g — go(b), and g = g c (b ), as 
is shown in Fig. 1. The four regions are as follows. 

A. When either 

0 < b < b 0 , g < g 0 (b) , (15) 

or 

b > b 0 , g < g c [b) , (16) 


there exists a single fixed point {£ 3 , 23 } that is a stable focus, which degenerates to a stable node 
in the vicinity of the line go(b). 

B. When either 

0 < b < - , g 0 (b) < g < 0 , (17) 

e 

or 

- <b<b 0 , g 0 (b) < g < g c (b) , (18) 

e 

there are three fixed points, an unstable focus {xi,z*}, a saddle and a stable node 

{ x 3i ^ 3 }- 

C. For 

b > - , g c (b ) < g < 0 , (19) 

e 


there is an unstable focus {x\, z *} and there appears a limit cycle around this focus. 

D. When 

0 < 6 < - , 0 <g< g c (b ) , 

e 

there are two fixed points, a saddle {x* 2 , z 2 } and a stable node {£ 3 , Z 3 }. 

E. Finally, when either 

0 < b < - , g> g c (b) , 


or 


b> ~, g >0 , 
e 


( 20 ) 


( 21 ) 

( 22 ) 


there are no fixed points. 

The transformation of fixed points, when moving from one region to another, can be illustrated 
by the corresponding bifurcation paths. 


A —> B —» D —» E: this bifurcation path is shown in Fig. 2, when moving from region 
A, through B and D, to region E along the line b = 0.2, for which 0 < b < 1/e. Then 
g 0 (b = 0.2) = —0.0194 and g c (b = 0.2) = 0.276. When the stable focus {£ 3 , £ 3 } in the region A 
approaches go(b), it degenerates to a stable node at g n = —0.0471 and then survives in the region 
B where, in addition, there appear an unstable focus {xl,z^} and a saddle Moving 

further from the region B to the region D , there remain only the stable node {£ 3 ,^ 3 } and the 
saddle {£ 2 , z* 2 }. Going from region D to region E by crossing g c {b), all fixed points cease to exist. 
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A — y B —>■ C —> E: this bifurcation path along the line b = 0.4, such that 1/e < b < b 0 = 
0.47, is illustrated in Fig. 3. Here g 0 (b = 0.4) = —0.0552 and g c (b = 0.4) = —0.0294. Moving 
upward in the region A, the stable focus {^ 3 ,^ 3 } becomes a stable node at g n = —0.0849. The 
stable node {£ 3 , z^\ survives into region B, where there appear two more fixed points, an unstable 
focus {x\,zf} and a saddle {x 2 ,z 2 }. In region C, the unstable focus {x\,z\} remains and there 
appears a limit cycle around it. In the region E, where g > 0, there are no fixed points. 

A — y C — y E: the bifurcation path along the line b = 1 > b 0 = 0.47 is given in Fig. 4. 
Then g c {b = 1) = —0.1769. The stable focus {x^,z^} in region A transforms in region C into 
an unstable focus {xl,zl} and there appears a limit cycle. The characteristic exponents at the 
critical line g c (b) are purely imaginary. No fixed points exist in the region E. 

A — y B — y C: the bifurcation path along the line g = —0.03 is presented in Fig. 5. The 
critical lines are crossed at the points 61 = b 0 (g = —0.03) = g^ig = —0.03) = 0.2718 and 
b 2 = b c (g = —0.03) = g~ l {g = —0.03) = 0.4007. The stable focus {^ 3 , 23 } in region A becomes 
a stable node at point b n = 0.1242 that survives into region B , where an unstable focus {x\,z\} 
and a saddle { x £, z%} appear. At the boundary b 2 , the stable node {X 3 , z 3 } and a saddle { xz%} 
merge together, while the unstable focus {x\,zf} continues to the region C, where a limit cycle 
arises. 

A — y C: the bifurcation path along the line g = —0.2 is shown in Fig. 6. The critical line 
g c (b) is crossed at the point b 2 = b c (g = —0.2) = g~ x (g = —0.2) = 1.1864. Here, the stable focus 
{£ 3 , 23 } in region A transforms into an unstable focus in region C, where also a limit 

cycle appears. 

The limit cycle, representing periodically occurring bubbles and crashes in asset prices, ap¬ 
pears in two cases. One is the path A — y C, with crossing the critical line g c (b ), accompanied by 
the bifurcation stable focus ==>- unstable focus + limit cycle, which is a typical Hopf bifurcation. 

The second case happens when crossing the line g c (b ) in the path B — y C, with the bifurca¬ 
tion stable node + saddle + unstable focus ==>- unstable focus + limit cycle. 

4 Periodically collapsing bubbles 

4.1 Different types of bubbles 

We find periodically sharply increasing and fast decreasing prices in region C. The amplitude, 
periodicity, as well as the shape of the price trajectories depend on the parameters b and g. 
Figure 7 shows the behavior of the asset price x(t) and bond prices z(t), where x(t) and z(t) 
are out of phase: an increase in the asset price is accompanied by a decrease in the bond price 
and vice-versa. In this parameter regime, the price rise and decay are symmetric, representing a 
typical business cycle regime. In contrast, Fig. 8 shows a parametrization, for which the price 
of the asset is exhibiting a bubble-like trajectory followed by a sharp faster correction. The 
prices remain periodic and the structure is that of periodically collapsing bubbles. We find that 
the asymmetry increases when g approaches the line g — 0, as is demonstrated in Fig. 9. The 
distance between the bubbles increases when approaching the line g c (b) with 1/e < b < b 0 , as 
illustrated in Figs. 10 and 11. 

We find that the maximum of the asset price x(t max ) occurs at a time t max that always precedes 
the time tf mn , where the bond price z(t z rnin ) has its minimum. The time lag At = t z min — t max 
depends on the parameters b and g. Rather than considering the absolute value of At, we 
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quantify its relative value reduced to the bubble width w measured at the half amplitude of the 


corresponding bubble. Varying the parameters b and g inside the region C, we find that the 
relative asset-bond time lag (A t/w) x 100% is of order 10%. 

An economic interpretation is as follows. The price rise during a bubble reflects an exuberant 
market that anticipates a booming economy. Two mechanisms concur to drive the bond prices 
lower, and thus the interest rate higher. As the economy is booming, it needs investment for 
growth. There is thus a scarcity of money, which makes its cost larger. In other words, the 
cost of credit increases and the interest rates rise. At the macro level, the central bank is also 
often tempted to moderate the exuberance of the bubble by making more expensive the access 
to credit, thus increasing the short-term interest rate that it can control. The fact, that the 
maximum of the asset price x(t max ) occurs at a time t max that always precedes the time t z min , 
where the bond price ^ ias its minimum, means that the credit dynamics and/or central 

bank intervention lags behind the price dynamics. This is in agreement with detailed analyses 
that have shown that stock markets in the US in general lead the Fed rate interventions |43 3 41J- 
If we consider that the dotcom bubble developed from Jan. 1997 to March 2000, covering roughly 
39 months, the relative lag of 10% we find in our model translates into about 4 months calendar 
time, which is in agreement with the estimation of the lag of about a quarter [431144] . 

4.2 Least stable phase during bubbles 

As seen in Figs. 8-11, the rising part of a bubble phase is characterized by a fast acceleration, 
which will be quantified precisely below. Given this strong growth regime, one can expect the 
prices to be more and more susceptible to external influences and shocks. Does the maximum 
susceptibility occur during the bubble ascent, at its maximum, or later during the crash? 

In the case of a strong asymmetry between the bubble ascent phase and the sharp crash 
correction, it is tempting to use an analogy with phase transitions and critical phenomena, where 
the point of phase transitions would correspond to the asset price peak. Then, one could expect 
that the response functions [T5j would be maximal at the points of phase transitions, where the 
system is supposed to be unstable [Ml31]. But our dynamical system is globally stable as a 
result of the feedback loops between the asset and bond prices. This translates into a globally 
stable periodic attractor, notwithstanding the very nonlinear nature of the oscillations, which 
we have termed “periodically collapsing bubbles” on reference to their shapes. Thus, the correct 
way to quantify the susceptibility is via a local stability measure called the relative expansion 
exponent [481149] . defined as the sum of the local Lyapunov exponents: the larger the expansion 
exponent, the lower the system stability. 

For the dynamical system (HUD , the local expansion exponent is defined as 118.TO' 



(23) 


where J(x,z ) is the Jacobian matrix 



1 dfi dfi 

dx dz 


(24) 


dh d/2 

dx dz 


for the dynamical system (fTOl) with 




Numerical investigation shows that the expansion exponent A is always negative, which confirms 
that the dynamical system is always stable in the sense that the dynamics remains bounded 
and non-chaotic. The largest value of A corresponds to the least stable state. It turns out that 
the market, described by the system of Eqs. m, is the least stable not where the asset price 
reaches its maximal value. This is illustrated in Fig. 12 for the same parameters (b = 0.4 and 
g = —0.029) used in Fig. 11. The asset price reaches its first maximum x(t max ) = 61.717 at 
tmax = 122.89, where the expansion exponent presents a local minimum A (t max ) = —0.878. The 
algebraic maximum of the expansion exponent A(t\) = —0.864 occurs at the time t\ = 121.48, 
which is slightly earlier than t max . Thus, the asset price is the least stable not at its maximum 
but just before it. 

4.3 Quantification of exploding bubbles 

We have found that the asset price x(t) can be well approximated during the bubble explosion 
by the following analytical form 

x ‘ pp{t) = exp {(iyvjy} ’ (25) 

where t\ is the point of maximum of the expansion exponent A(t). 

Our numerical investigations lead us to conclude that the exponents a and (3, for the values 
of b and g in the left corner of the region (7, between region B and E, where 1/e < b < b 0 and 
g > g c (b), are close to 

a = l , P = \. (26) 

The numerical coefficients c± and C 2 are positive parameters that depend on the specific values of 
b and g. For instance, for b = 0.4 and g = —0.029, we have ci = 2.8, C 2 = 1.5, and t\ = 121.48. 

It is interesting to find out whether there exists a critical region where the parameters a 
and (3 are invariant, being close to a = 0.4 and = 0.2, respectively. We have accomplished a 
detailed analysis over the whole region C. It turned out that approximation (|25j) describes very 
well the transitional behavior of the beginning of the bubble bursting for the parameters b and 
g in the vicinity of the double critical point {b = 1/e , g = 0}, which is the intersection point of 
the critical lines b = b c (g) and g — 0, where the exponents a and f3 play the role of a kind of 
critical indices. 

Moreover, the super-exponential form of the initial stage of bursting bubbles is well repre¬ 
sented by approximation (j25j) in a finite region of the parameters b and g corresponding to the 
triangle, where 1/e < b < b 0 and g > g c (b). In this critical region, the exponents a and [3 are in¬ 
variant, and only the parameters C\ and c 2 have to be adjusted. For example, let us take b = 0.38 
and g = —0.0117, for which g c = —0.01174. Then the asset price reaches its first maximum 
x(tmax) = 393.5 at tmax ~ 501.5, where the expansion exponent possesses its local minimum 
A (tmax) = —0.9607, while the absolute maximum of the expansion exponent A{t\) = —0.9535 
occurs at the time t\ ~ 499.74 < t max . Then the critical exponents a = 0.4 and /3 = 0.2 can be 
kept invariant, and only the parameters c\ ~ 1.5 and C 2 ~ 2.8 are varied. 

In order to give a deeper understanding justifying why and where the super-exponential form 
Q25J) serves as a good description of the transition from a smooth, almost unchangeable, behavior 
of the asset price to the start of a bursting bubble, let us study the ratio R = width/period 
defined as the ratio of the temporal bubble width, measured at half of its amplitude, to the 
temporal period between the bubbles. It is reasonable to expect that, for genuine bubbles, the 
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ratio R has to be small, representing sharp bubble bursts and crashes interspersed within a 
smooth market behavior. Our numerical investigation shows that the super-exponential form 
([25]) . with invariant exponents a = 0.4 and /3 = 0.2, is valid as long as R -C 1, which occurs in 
the triangle region where 1/e < b < bo ~ 0.47 and g c < g < 0. 

The super-exponential behavior (j25l) is reminiscent of the many reported empirical cases of 
financial bubbles [311501 - 154] . whose structure can serve as a precursor of the appearing bubble [28]. 
We note that the transient super-exponential explosive bubble trajectory is here stronger than the 
simple power law singularities discussed elsewhere [521 - 154] as form (125|) represents an intermediate 
asymptotic with an essential singularity behavior. Figure 13 demonstrates the excellent accuracy 
of the approximation (]25l) . as compared to the accurate “exact” numerical solution x{t ) in the 
region of its accelerated growth. 


4.4 Period and amplitude of the periodically collapsing bubbles 

The shape of the bubbles and the distance between them strongly depends on the closeness of 
the parameters to the critical lines g = g c (b) and g = 0. This results from the existence of a 
stable node in the region B, which influences the length of the plateaus between the bubbles in 
the region C , as is seen in Figs. 8 and 11. Increasing the distance from the boundary, separating 
the regions B and C, that is, increasing b with fixed g , shortens the plateau length between two 
successive bubbles, but does not influence much the bubble amplitude. This is illustrated in Fig. 

14. The level of the plateaus does not depend on initial conditions, as is demonstrated in Fig. 

15. 

The boundary g c (b ) plays the role of a critical line at which the plateau length diverges, 
similarly to the divergence of the correlation length of statistical physics systems at a phase 
transition line j46[|47]. To describe this divergence, let us define the duration 

L = t n+ \ — t n (27) 


between two successive maxima of the asset price x{t n ) occurring at the times t n and t n+ \. In 
other words, L is the time interval between two successive bubbles. This interval, for a fixed g , 
essentially depends on the deviation 

A = b-b c {g) (28) 

of b from the boundary point b c (g) = g^^g). The critical behavior can be characterized by the 
power law 

L = L( A) oc A" 1 ' (A ->• +0) . (29) 

Hence, the critical index is given by the limit 


v = — lim 

A—H-0 


cl\wL 
din |A| 


(30) 


Figure 16a shows the behavior of lnL as a function of A, from where we determine v = 0.5. 

The other critical line is g = 0 at which the bubble amplitudes diverge. Above this line, 
the behavior of the asset price is exponentially inflating. Approaching the critical line g — 0 
from below, the bubble amplitude A = x(t max ) diverges, analogously to the behavior of the 
compressibility in statistical physics systems close to a phase transition line, as 


A = A(g) oc g 7 (g-> 0 ) . 


(31) 
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This defines the critical index 


7 = — lim 


d in A 


(32) 


o d ill | <71 

The dependence of in A versus in |c/| is shown in Fig. 16b, from which we determine 7 = 1. 
Table 1 illustrates the fast divergence of the bubble amplitude A, when decreasing g with fixed 
6 = 1, the number N of bubbles in the time interval [0,100], and the bubble width w measured 
at half-amplitude. The bubble width, as g —> —0, saturates to 0.7, while its amplitude diverges. 


5 Concluding remarks 

We have identified a specific structure of feedback loops between asset prices and bond prices 
that lead to the regimes of periodically explosive bubbles followed by crashes. While asset and 
bond prices tend to mean revert to their respective fundamental prices, we have argued that 
the later are not fixed but depend themselves on the state of the economy represented by the 
instantaneous asset and bond prices. This results in a rich nonlinear coupled dynamics of asset 
and bond prices. The two key parameters b and g represent respectively the sensitivity of the 
fundamental asset price on past asset and bond prices and of the fundamental bond price on 
past asset prices. They thus capture the competition of the allocation of capital between assets 
and bonds that develops as the two exhibit diverging trajectories. We stress that the regime of 
periodically collapsing bubbles does not require heterogeneous agents or sudden changes in market 
characteristics. Periodically occurring bubbles and crashes can be treated as natural phenomena 
objectively arising in stock markets due to the delay of the impact of present valuations on future 
fundamental prices and the nonlinear couplings. 

The coupled nonlinear equations (ITU]) have been motivated as the description of the dynamics 
of a representative asset price coupled with a bond price. There are other possible interpretations. 
For instance, we can think of the variable x(t) as the monetary base of an economy and z(t) as 
the spread (or difference) between the yield (interest rate) of bonds with long term maturity (say 
10 years) and the interest rate of bonds with short term maturity (say 3 months). The spread 
z(t) provides a measure of the market expectation on the growth potential of the economy. The 
larger z(t), the larger the expected growth rate. When x is increasing and z is large, this means 
that short-term interest rates are low (cheap funding and borrowing) and long term interest 
rates are increasing (expectation of a growing economy with growing returns). This boosts the 
monetary base on the anticipation of a growing economy and as a result of the growth of the 
economy that needs financing. This is the mechanism by which a larger z promotes a future 
larger monetary base x. The multiplicative nature of the equations (fTOj) capture the fact that 
the positive effect of z on x is all the larger, the larger are x and z. But as the monetary base 
x increases, it generally grows faster than the real economy, leading to misallocation of funding 
to non-performing industry sectors. This then tends to push z down, on the expectation that 
long-term interest rates will become smaller as a result of the slowing down of the economy. 
Also, the economy in its boom phase is over heating and usually requires the intervention of the 
monetary authorities, who are led to increase the short term interest rates (the cost of short¬ 
term borrowing) to stabilize the economy. The larger x, the more negative its impact on z, 
hence on its limiting fundamental value. This is captured by the dynamic equations (TTTTj) . Thus, 
the regimes of parameters in region C, that we have reported in this article, demonstrate the 
difficulties in stabilizing an economy with the usual monetary tools consisting in attempt to 
control the monetary base and the interest rates. Due to intrinsic delay and nonlinear feedbacks, 
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such attempts can lead to unintended cycles of booms and busts, as documented for the period 
from 1980 to present and its role in the great crisis of 2008, following great recession and on-going 
global crises [55] . 
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Figure captions 


Figure 1. Region of existence of stationary solutions corresponding to the fixed points of 
the evolution equations (TTOj) for the asset price x and bond price z. 

Figure 2. Bifurcation path A —» B —» D —> E. Here b = 0.2 < 1/e is fixed and g 
varies: (a) asset price x*; (b) bond price z*. The fixed point {xg, z^} (solid line) is a stable 
focus transforming into a stable node (solid line with dots) at g n ~ —0.0471. The point exists 
for — oo < g < g c . The fixed point {x 2 , z 2 \ (dashed-dotted line) is a saddle, which exists for 
go < g < g c , and the fixed point {x\,z\} (dashed line) is an unstable focus, which exists for 
g 0 < g < 0. At g — g 0 , the fixed points {x^, z^} and {x^, z 2 } coincide. At g — g c , the fixed points 
{x?j, } and {x 2 ,z 2 } coincide. When g —> —0, then x\ —> +00 and z* —> 0. When g —> —00, 

then X3 —* 1 and z\ —> 0. When g > g c , fixed points do not exist. 

Figure 3. Bifurcation path A —>■ B —>■ C —>■ E. Here 1/e < b = 0.4 < b 0 and g varies: (a) 
asset price x*; (b) bond price z*. The fixed point {x^, Zg} (solid line) is a stable focus transforming 
to a stable node (solid line with dots) at g n ~ —0.0849. The point exists for —00 < g < g c < 0. 
The fixed point {x 2l z 2 } (dashed-dotted line) is a saddle, which exists for g 0 < g < g c , and the 
fixed point {x^,;^} (dashed line) is an unstable focus, which exists for g 0 < g < 0. At g = g 0 , 
the fixed points {xg,^} and {x 2 ,z 2 } coincide. At g = g c , the fixed points {x^Zg} and {x 2) z 2 } 
coincide. When g —> —0, then x\ —* +00 and z* 0. When g —>- — 00 , then X3 —> 1 and Z3 —> 0. 
When g > 0, fixed points do not exist. 

Figure 4. Bifurcation path A —> C —)• E. Here b = 1 > bo is fixed and g varies: (a) 
asset price x*; (b) bond price z*. The fixed point {xg,Xg} (solid line) is a stable focus, and the 
fixed point {xg,^} (dashed line) is an unstable focus. At g = g c ~ —0.1769, the stable focus 
transforms into an unstable focus with Xg = x\ = e 2 , z% = z 2 = 2 /(be 2 ), the Lyapunov exponents 
being = 0. 

Figure 5. Bifurcation path A —> B —)■ C. Here g = —0.03 is fixed and b varies: (a) asset 
price x*; (b) bond price z*. The fixed point {xg,£g} (solid line) is a stable focus transforming 
to a stable node (solid line with dots) at b n zz 0.1242. The point exists for 0 < b < &2, where 
b 2 ~ 0.4007. The fixed point {x 2) z 2 } (dashed-dotted line) is a saddle, which exists for b\ < b < b 2: 
where b\ zz 0.2718. The fixed point {xl,zl} (dashed line) is an unstable focus, which exists for 
b > bi. At b = bi, the fixed points {x\,zl} and {x^,^} coincide. At b = b 2j the fixed points 
{xg,Zg} and {xg,^} coincide. When b = 0, then Xg = 1 and Xg = e 9 . When b —> + 00 , then 
x* —> Too and z^ —> 0. 

Figure 6. Bifurcation path A —> C. Variation of fixed points under fixed g = —0.2 and 
changing b: (a) asset price x*; (b) bond price z*. The fixed point {xg,£g} (solid line) is a stable 
focus that at b = b 2 ze 1.1864 transforms into an unstable focus {x^, z *} (dashed line). At b = b 2l 
the fixed points {xg,£g} and {x^,Zg} coincide, with the Lyapunov exponents being 9fJAi i2 = 0. 

Figure 7. Behavior of the solutions for the asset price and bond price for the parameters 
b = 0.5, g = —0.083, and the initial conditions xq = 3 and zq = 0.1. Here, g c ~ —0.083056. 
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(a) Asset price x(t) for t G [0,600]; (b) bond price z(t) for t G [0,600]; (c) a single zoomed out 
bubble of the asset price x(t) for t G [495,510]; (d) a zoomed out negative bubble of the bond 
price z(t) for t G [495,510]. 

Figure 8. Asset price and bond price for the parameters 6 = 0.4, g = —0.0294, and the initial 
conditions x 0 = 3 and z 0 = 0.1. Here, g c « —0.029424. (a) Asset price x(t) for t G [0,1400]; (b) 
bond price z(t) for t G [0,1400]; (c) a zoomed out bubble of the asset price x(t) for t G [495, 510]; 
(d) a zoomed out negative bubble of the bond price z(t) for t G [495, 510]. 

Figure 9. Change in the behavior of the asset price x(t), with the initial conditions Xq = 3, 
z 0 = 0.1, parameters b — 1, and varying g: (a) g = —0.176 > g c , where g c (b) = —0.176862964. 
The asset price x(t) oscillates without convergence as t —» oo; (b) a bubble of x(t) for t G [69, 75], 
with the same parameters as in (a); (c) g = —0.05. The asset price x(t) displays periodic bubbles 
and crashes, as t —* oo. (d) a bubble of the asset price x(t) for t G [69,75], under the same 
parameters as in (c). 

Figure 10. Logarithmic behavior of the asset price for the parameters b — 1, g — —0.001, and 
the initial conditions Xq = 3 and Zq = 0.1. (a) The function lnx(t) oscillates without convergence 
for t —» oo; (b) a bubble of In x(t) for t G [ 66 , 84] under the same parameters as in (a). 

Figure 11. Logarithmic behavior of the asset price x(t) for the parameters b = 0.4, g = 
—0.029, and the initial conditions x 0 = 1 and z 0 = 0.1. Here, g c (b) = —0.029424. (a) The 
function In x(t) oscillates without convergence as t —* oo; (b) a bubble of hax(t) for t G [0,160], 
with the same parameters as in (a). 

Figure 12. Expansion exponent A (t), with the parameters b = 0.4, g = —0.029, and the 
initial conditions x 0 — 1 and z 0 = 0.1, for different temporal scales: (a) t G [100,1000]; (b) 
t G [110,130]. At t max ~ 122.89, the asset price x(t) has its first local maximum x max = 61.7171, 
which corresponds to the first local minimum of A (t). 

Figure 13. Comparison of the numerical solution for the asset price x(t) (solid line) and 
its approximation x app (t ) (dashed-dotted line), with initial conditions xq — 1 and zo = 0 . 1 , for 
different parameters b and g, where 1/e < b < b 0 « 0.47, and for different intervals of time: 
(a) b = 0.4, g = —0.029 > g c — —0.02943, and t G [0,121]; (b) b and g , as in Fig. 13a, but 
t G [110,121]; (c) b = 0.38, g = —0.0117 > g c = —0.01174, and t G [0,499.6]; (d) b and g , as in 
Fig. 13c, but t G [489,499.6]; 

Figure 14. Change in the behavior of the asset price x(t), with fixed g = —0.03, when 
varying b in the vicinity of the boundary point 62 = 0.400691 separating the regions with different 
numbers of fixed points. The initial conditions are x 0 — 1 and z {) — 0 . 1 . (a) Asset price x(t) for 
b = 0.401 > & 2 (solid line) and for b = 0.41 (dashed-dotted line) for t G [0,200]; (b) asset price 
x(t) for b = 0.401 > 62 (solid line) and for b = 0.501 (dashed-dotted line) for t G [0,200]; (c) 
asset price x(t) for 6=1 (solid line), for 6 = 2 (dashed line), and for 6 = 10 (dashed-dotted line), 
with t G [0, 20]. 

Figure 15. Asset price x(t) for the parameter g = —0.03, but with different 6 and for different 
initial conditions {xo,^o}- ( a ) Asset price x(t) for 6 = 0.4007 > 6 2 (solid line), 6 = 0.4006 < 6 2 
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(dashed-dotted line), and the initial conditions x 0 — 1 and z 0 = 0.1. Here, 6 2 = 0.400691 is 
the boundary point, such that for b > 6 2 , there is an unstable focus and a limit cycle, while 
for b < b- 2 , there are three fixed points. The fixed point { x 3 = 2.928, Zi = 0.916} is a stable 
node, with the Lyapunov exponents {Ai = —0.903, A 2 = —0.022}. The fixed point {x 2 = 
3.08, Z 2 = 0.912} is a saddle, with the Lyapunov exponents {Ai = —0.899, A 2 = +0.022}, and 
the fixed point {x\ = 58.26, z± = 0.174} is an unstable focus, with the Lyapunov exponents 
{Ai = A2 = +1.03 — il.72}. (b) Asset price x(t ) for the parameters g = —0.03, b = 0.4007 > 6 2 , 
and different initial conditions: {x 0 — l,z 0 — 0.1} (solid line), {xq = 5, z 0 = 0.5} (dashed line), 
and {xo = 0.01, z 0 = 50} (dashed-dotted line). 

Figure 16. Critical behavior of bubble characteristics. The temporal interval between two 
neighboring bubbles is L. The bubble amplitude is A. The distance from the critical line is 
A = b - b 2 in (a) and g in (b). The initial conditions are taken as Xq = 1 and Zo = 0.1. (a) 
Interval L in logarithmic units, as a function of log-distance ln|A| = In |6 — b 2 \ to the critical 
line. Here, g = —0.03 and & 2 = 0.400691. (b) Bubble amplitude A in logarithmic units, as a 
function of ln|g| measuring the distance from the critical line g = —0, with fixed 6 = 1. Here 
g c = —0.1769, so that \i\g c = —1.7322. 
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Table Caption 


Table 1. Fast divergence of the bubble amplitude A, when decreasing g, with fixed 6 = 1. 
N is the number of bubbles in the time interval [0,100], and w is the bubble width measured at 
half-amplitude. 
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Figure 1: Region of existence of stationary solutions corresponding to the fixed points of the 
evolution equations for the asset price x and bond price z. 
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Figure 2: Bifurcation path A — y B — y D — y E. Here b = 0.2 < 1/e is fixed and g varies: (a) 
asset price x*; (b) bond price z*. The fixed point {xg, 2:3} (solid line) is a stable focus transforming 
into a stable node (solid line with dots) at g n ~ —0.0471. The point exists for — 00 < g < g c . 
The fixed point {x^, A/} (dashed-dotted line) is a saddle, which exists for g 0 < g < g c , and the 
fixed point {xl,z^} (dashed line) is an unstable focus, which exists for g 0 < g < 0. At g — g 0 , 
the fixed points |x*, z*} and {x^,z%} coincide. At g — g c , the fixed points {x^,z^} and {x%,z%} 
coincide. When g —y —0, then x* —y +cxd and z\ —y 0. When g —y — 00 , then X3 —y 1 and z% —y 0. 
When g > g c , fixed points do not exist. 
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Figure 3: Bifurcation path A —y B —y C —y E. Here 1/e < b = 0.4 < bo and g varies: (a) 
asset price x *; (b) bond price z*. The fixed point {£3, Z3} (solid line) is a stable focus transforming 
to a stable node (solid line with dots) at g n ~ —0.0849. The point exists for —00 < g < g c < 0. 
The fixed point {x%,z%} (dashed-dotted line) is a saddle, which exists for g 0 < g < g c , and the 
fixed point {£3,23} (dashed line) is an unstable focus, which exists for g 0 < g < 0. At g — g 0 , 
the fixed points {re) 1 , z*} and {x^z^} coincide. At g = g c , the fixed points {0:3, ^3} and {x^z^} 
coincide. When g —y —0, then x\ —y +cxd and z* —y 0. When g —y — 00 , then x^ —y 1 and £3 —y 0. 
When g > 0, fixed points do not exist. 
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Figure 4: Bifurcation path A —y C —> E. Here b = 1 > &o is fixed and g varies: (a) asset price 
x *; (b) bond price z*. The fixed point {x^,z^} (solid line) is a stable focus, and the fixed point 
{x\,z \(} (dashed line) is an unstable focus. At g = g c ~ —0.1769, the stable focus transforms 
into an unstable focus with = x\ = e 2 , ^3 = z* = 2 /(be 2 ), the Lyapunov exponents being 

KAi,2 = 0. 
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Figure 5: Bifurcation path A — y B — y C. Here g = —0.03 is fixed and b varies: (a) asset price 
x*; (b) bond price z*. The fixed point {xg, £3} (solid line) is a stable focus transforming to a stable 
node (solid line with dots) at b n w 0.1242. The point exists for 0 < b < 62, where 62 ~ 0.4007. 
The fixed point {x^z^} (dashed-dotted line) is a saddle, which exists for b\ < b < fe 2 , where 
61 ~ 0.2718. The fixed point {x^,;?*} (dashed line) is an unstable focus, which exists for b > b\. 
At b = 61, the fixed points {xj, z*} and { x 2, ^2} coincide. At b = & 2 , the fixed points {xg, Zg} and 
{x^, ^2} coincide. When b = 0, then Xg = 1 and Zg = e 9 . When b —y +cx), then x^ —y +cx) and 
z{ -> 0. 
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Figure 6: Bifurcation path A —» C. Variation of fixed points under fixed g = —0.2 and changing 
b\ (a) asset price x *; (b) bond price z*. The fixed point {£3, £3} (solid line) is a stable focus that 
at b = b -2 ~ 1.1864 transforms into an unstable focus {£3, z 3} (dashed line). At b = b 2 , the hxed 
points {^3,^3} and {£3,2*} coincide, with the Lyapunov exponents being 9ftAi j2 = 0. 


25 





Figure 7: Behavior of the solutions for the asset price and bond price for the parameters b = 0.5, 
g = —0.083, and the initial conditions xq = 3 and zq = 0.1. Here, g c ~ —0.083056. (a) Asset 
price x(t) for t G [0,600]; (b) bond price z(t ) for t G [0,600]; (c) a single zoomed out bubble of 
the asset price x(t) for t G [495,510]; (d) a zoomed out negative bubble of the bond price z(t) 
for t G [495, 510]. 
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Figure 8: Asset price and bond price for the parameters b = 0.4, g = —0.0294, and the initial 
conditions x Q = 3 and z 0 = 0.1. Here, g c ~ —0.029424. (a) Asset price x(t) for t G [0,1400]; (b) 
bond price z(t ) for t G [0,1400]; (c) a zoomed out bubble of the asset price x(t) for t G [495, 510]; 
(d) a zoomed out negative bubble of the bond price z{t ) for t G [495, 510]. 


27 























Figure 9: Change in the behavior of the asset price x(t), with the initial conditions Xq = 3, 
z 0 = 0.1, parameters 6 = 1, and varying g: (a) g = —0.176 > g c , where g c (b) = —0.176862964. 
The asset price x(t) oscillates without convergence as t —> oo; (b) a bubble of x(t) for t G [69, 75], 
with the same parameters as in (a); (c) g = —0.05. The asset price x(t) displays periodic bubbles 
and crashes, as t —> oo. (d) a bubble of the asset price x(t) for t G [69,75], under the same 
parameters as in (c). 
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Figure 10: Logarithmic behavior of the asset price for the parameters b — 1, g — —0.001, and the 
initial conditions Xo = 3 and Zq = 0.1. (a) The function lnx(f) oscillates without convergence for 
t —> oo; (b) a bubble of lnx(f) for t 6 [66, 84] under the same parameters as in (a). 
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Figure 11: Logarithmic behavior of the asset price x(t) for the parameters b = 0.4, g = —0.029, 
and the initial conditions Xq = 1 and Zq = 0.1. Here, g c (b) = —0.029424. (a) The function lnx(t) 
oscillates without convergence as t —> oo; (b) a bubble of In x(t) for t e [0,160], with the same 
parameters as in (a). 
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Figure 12: Expansion exponent A (t), with the parameters b = 0.4, g = —0.029, and the initial 
conditions x 0 = 1 and z 0 = 0.1, for different temporal scales: (a) t £ [100,1000]; (b) t £ [110,130]. 
At t max ~ 122.89, the asset price x{t ) has its first local maximum x max = 61.7171, which 
corresponds to the first local minimum of A(f). 
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Figure 13: Comparison of the numerical solution for the asset price x(t) (solid line) and its 
approximation x app (t ) (dashed-dotted line), with initial conditions x 0 — 1 and z 0 = 0.1, for 
different parameters b and g, where l/e<6<6 0 ~ 0.47, and for different intervals of time: 
(a) b = 0.4, g = —0.029 > g c — —0.02943, and t G [0,121]; (b) b and g, as in Fig. 13a, but 
t G [110,121]; (c) b = 0.38, g = —0.0117 > g c = —0.01174, and t G [0,499.6]; (d) b and g, as in 
Fig. 13c, but t G [489,499.6]; 
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Figure 14: Change in the behavior of the asset price x(t), with fixed g = —0.03, when varying 
6 in the vicinity of the boundary point b 2 = 0.400691 separating the regions with different 
numbers of fixed points. The initial conditions are x Q — 1 and Zo = 0.1. (a) Asset price x{t) for 
b = 0.401 > b 2 (solid line) and for b = 0.41 (dashed-dotted line) for t G [0,200]; (b) asset price 
x(t) for b = 0.401 > b 2 (solid line) and for b = 0.501 (dashed-dotted line) for t G [0,200]; (c) 
asset price x{t) for 6 = 1 (solid line), for 6 = 2 (dashed line), and for 6 = 10 (dashed-dotted line), 
with t G [0, 20]. 
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Figure 15: Asset price x(t) for the parameter g = —0.03, but with different b and for different 
initial conditions {xo, Zo}- (a) Asset price x(t) for b = 0.4007 > b 2 (solid line), b = 0.4006 < b- 2 
(dashed-dotted line), and the initial conditions xq — 1 and zq = 0.1. Here, b 2 = 0.400691 is 
the boundary point, such that for b > b 2 , there is an unstable focus and a limit cycle, while 
for b < b 2 , there are three fixed points. The fixed point {£3 = 2.928, Z\ = 0.916} is a stable 
node, with the Lyapunov exponents {Ai = —0.903, X 2 = —0.022}. The fixed point {x 2 = 
3.08, z 2 = 0.912} is a saddle, with the Lyapunov exponents (Ai = —0.899, A 2 = +0.022}, and 
the fixed point {x\ = 58.26, z\ = 0.174} is an unstable focus, with the Lyapunov exponents 
{Ai = A 2 = +1.03 — 71.72}. (b) Asset price x(t) for the parameters g = —0.03, b = 0.4007 > b 2 , 
and different initial conditions: {^o = l,z 0 — 0.1} (solid line), (x 0 = 5 ,z 0 — 0.5} (dashed line), 
and {xq = 0.01, z 0 = 50} (dashed-dotted line). 
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Figure 16: Critical behavior of bubble characteristics. The temporal interval between two neigh¬ 
boring bubbles is L. The bubble amplitude is A. The distance from the critical line is A = b — b 2 
in (a) and g in (b). The initial conditions are taken as xq = 1 and zq — 0.1. (a) Interval L 
in logarithmic units, as a function of log-distance In |A| = In \ b — £> 2 1 to the critical line. Here, 
g = —0.03 and 62 = 0.400691. (b) Bubble amplitude A in logarithmic units, as a function of 
In | g | measuring the distance from the critical line g = —0, with fixed b = 1. Here g c = —0.1769, 
so that In g c = —1.7322. 
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